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Abstract 

Recently proposed by Battistel et al. "predictive formulation of the NJL model" is dis- 
cussed and its connection with the differential regularization is noted. The principal problem 
of this formulation is a non-physical singularity (Landau pole) in meson propagators. A modi- 
fication of the formulation, which is free of the Landau pole and conserves main features of the 
approach, is proposed. 
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Effective models in the non-perturbative region are a considerable part of the modern strong- 
interaction theory. One of the most successful effective models of strong interaction of light hadrons 
is the Nambu-Jona-Lasinio (NJL) model with quark content. Since the NJl model in the leading 
approximation includes the quark loops and is based on non-renormalized four-fermion interaction, 
the essential aspect of this model is a regularization, which on the widely current opinion constitutes 
a part of the definition of the NJL model. Though predictions of the NJL model for commonly 
used regularizations (such as the cutoff in momentum space, the Pauli-Villars regularization or the 
proper-time regularization) are very similar (see, e.g., [2] for review), nevertheless such dependence 
on the regularization prescription cannot be satisfactory from the general theoretic point of view. 

A very interesting attempt was made recently in the work of Battistel et al [1] for releasing 
the NJL model from the regularization dependence. In this work a method for the definition of 
one-loop Green functions of the NJL model was proposed (see also foregoing works [3]). The idea of 
the method is to avoid the explicit evaluation of divergent integrals with any specific regularization. 
The finite parts are separated of the divergent ones and are integrated without regularization. Then 
the NJL model becomes predictive in the sense that its consequences do not depend on the specific 
regularization of divergent integrals. An important result of work [1] is a proof of the fulfilment 
of all symmetry constraints on leading-order Green functions. A parameter choice leads to the 
reasonable values of the constituent quark mass and other model parameters. 

In the present paper some features of the calculational scheme of work [1]^ are analyzed, and 
a connection of the BDK approach with the differential regularization is pointed. The principal 
problem of the BDK approach is the singularity of meson propagators in the Euclidean momentum 
region. The presence of this singularity (Landau pole) prevents mcson-loop calculations and, there- 
fore, makes impossible any calculations of corrections to the leading approximation. A modification 
of the scheme, which is free of Landau pole and conserves the main features of the approach of [1], 
is proposed. 

Following [1] we consider the simplest physically non-trivial variant of the NJL model with 
chiral symmetry SUvi2) x SUa{2). The model Lagrangian is 



Here ip is the quark field with ric colours, mo is the current quark mass, g is the coupling constant 
of dimension and T are Pauli matrices. 

The leading approximation of the model is the mean-field approximation, which coincides in the 
case with the leading order of 1/nc-expansion. All Green functions of the leading approximation 
are expressed in terms of quark one-loop integrals. A problem of calculations of these integrals is 
reduced to the definition of following five divergent integrals (see [1]): 



Here M is the dynamical (constituent) quark mass, which is a non-trivial solution of the gap 
equation for the NJL model. 

A basic device for the definition of integrals (2) and (3) in the BDK approach is an alge- 
braic identity for the propagator function (see equation (28) in work [1]). With the identity the 
divergent parts of integrals (2) and (3) have been expressed via three tcnsorial and two scalar 
external-momentum-indcpendent integrals which were treated in the sense of some unspecified reg- 
ularization. Then integrals (2) and (3) have been expressed in terms of these five integrals and two 
standard convergent integrals (see formulae of section III in work [1]). 

'^We shall name the calculational scheme of work [1] as the BDK approach, or the implicit regularization. 
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In connection with these expressions we draw attention to some general regularization-independent 
property of integrals (2) and (3). Namely, integrals {Ii; If} and {I2; I2',} are connected with fol- 
lowing relations 



and 



■.l2{ki,h; M^) (4) 

^^2 M'). (5) 

These relations are regularization independent and play a significant part in the approach. It is 
easy to verify the validity of these relations for traditional regularization schemes (four-dimensional 
cutoff, Pauli-Villars regularization, etc) by straightforward calculation, but their validity is not 
evident directly from above-mentioned formulae of work [1]. To prove these relations we should 
define derivatives of the external-momentum- independent integrals over M^. For this purpose we 
note that derivatives of logarithmically-divergent integrals are convergent integrals and can be 
calculated without any regularization. The calculations of these convergent integrals gives us zero 
value for derivatives of tensorial logarithmically-divergent integrals. Then it is easy to verify that 
the derivatives of quadratically-divergcnt integrals arc the corresponding logarithmically-divergent 
integrals. Taking into account these circumstances it is easy to prove equations (4) and (5) for 
BDK expressions of integrals (2) and (3). The significance of relations (4) and (5) are evident since 
they can be used for alternative definitions of quadratically divergent integrals Ii and /f without 
additional regularization. 

With obtained expressions for integrals (2) and (3) the authors of work [1] have analyzed the 
symmetry constraints on one-loop Green functions and have argued the necessity of consistency 
relations, which are concluded in zero value for tensorial external-momentum-independent integrals 
(see eq. (86) in work [1]). With these relations all symmetry properties of the theory (sTich as Furry 
theorem, Ward identities, etc) are fulfilled for all one-loop Green functions. This point is one of 
the important results of work [1]. From the point of view of above discussion the BDK consistency 
relations simply assert zero values of corresponding integration constants. 

The choice of parameters of the NJL model with Lagrangian (1) in the leading approximation 
is principally defined by two divergent integrals, namely, logarithmically-divergent integral I2 and 
quadratically-divergent integral Ii. Integral Ii is a part of the gap equation 

M = mo - SigUcM ■ h, (6) 

which defines dynamical (constituent) quark mass M. Integral I2 determines the structure of meson 
propagators. Both these integrals can easily be defined with differential regularization without using 
the above-mentioned algebraic identity for the propagator function. 

To define I2 one can use well-known trick, which goes back to Gelfand and Shilov [4] (see also 
[5]). Namely, let introduce new external variables 

p = ki-k2, P = ki + k2. (7) 

Derivatives dl2ldpfji, 812/ dP^ are convergent integrals, and their calculation gives us 

dh ^ iPu. , u{l - u) dh ^ . . 

dpj, 167r2yo M2-n(l-n)p2' dP^ 

Basing on these results and using identity p^-^ = \Pn^^i we naturally go to the following defini- 
tion: 

. / 9x i /"^ , , - m(1 - u)r?2 
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where Mq is the integration constant. This expression coincide with the BDK ones if constant Mg 
is related to the scalar external- momentum-independent integral Iiog (see eq. (40) in work [1]) as 

^^o. = -^log^ = /.(0). (10) 

An essential feature of the above calculation is the permutation of two limits - differentiation 
and regularization removing. Sure, such permutation is implied also in work [1] (in the calculations 
of finite parts). This permutation is an essence of Gelfand-Shilov differential regularization [4], 
which is based on the infinite differentiability of generalized functions. 

To define Ji we use regularization-independent relation (4), which gives (with taking into ac- 
count eq. (10)) 

A = ^(M^log§ + M^-Mf), (11) 

where Mf is the integration constant of eq. (4). By taking into account the consistency relations 
such definitions also corresponds to BDK ones. 

To calculate Mq and Mi and, therefore, to define Iq and Ii in full, we can use in the chiral limit 
(mo = 0) two regularization-independent relations of NJL model (see, e.g. [2]), namely 

t; = '-^m (12) 

where /tt is the pion-decay constant, and 

X = 2c3 = (13 
where x =< Ol^/S-^IO > and c is the quark condensate. We have from equation (12) 

'^(^'=16^[i^-/„ <flog(l-„(l-„)^)]. (14) 
Then, using equations (11) and (13), we obtain for M\ 

M^ = M^ + ^^-^. (15) 

^ Uc UcM ^ ' 

This equation (wc shall refer it as BDK equation) plays a key role in the BDK approach.^ 

At pion constant and quark condensate c being fixed equation (15) is the equation for 
dynamical quark mass M, where M\ is a parameter. In terms of new variable x = -(l^SzT^l^i ^ 
equation (15) can be written as 

- 3ax 2 = 0, (16) 

where a = ^""'^g^ii'* ' {M\ — '^^Jj )- Depending on the value of parameter a, three cases are possible: 

(1) at a < 1 equation (16) possesses one real negative root x\ < 0; 

(2) at a = 1 equation (16) possesses one negative root xi = —2 and one positive root X2 = 1. 

(3) at a > 1 equation (16) possesses one negative root xi < and two positive roots < a;2 < 1 
and xs > 1. 

It is clear, that at c < the second case (a = 1) only can be physically accepted. The authors 
of work [1] take only this situation as a solely possible choice for model parameters. The value of 
the dynamical quark mass, which corresponds to positive root X2 = 1 at ric = 3, is M = —1.873 c. 



An alternative way to introduce the implicit regularization can be based on an approach of work [6]. 
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At condensate value c = —250 MeV it gives M = 468 MeV. Coupling constant g can be defined 
from the well-known relation of the NJL model 



g = -M/2c^ (17) 

which is also regularization-independent (see, e.g., [2]). Noteworthy, the quark mass value and the 
coupling value depend on the quark condensate value only, and do not depend on the pion constant 
value. The last one defines the value of Mi, which is Mi = 879 MeV at fj^ = 93 MeV. 

It is interesting to compare the parameter values of the BDK approach with the parameter 
values of traditional regular izations. The value of quark mass (468 MeV) in the BDK approach 
is noticeably higher in comparison with the values in the four-dimensional momentum cutoff reg- 
ularization (236 MeV) and the Pauli-Villars regularization (240 MeV) (at given value of quark 
condensate -250 MeV). Apart from the quark mass dependence on the condensate value is quite 
different. In the BDK approach the quark mass is proportional to the condensate, whereas in 
4-momentum cutoff and Pauli-Villars regularization the quark mass increases at decreasing the 
absolute value of condensate. For instance, at condensate value c = —210 MeV the quark mass 
is M = 393 MeV in the BDK approach and M = 423 MeV in the foTir-dimcnsional momentum 
cutoff. At the same time, from the point of view of the phenomenology the parameters' values in 
BDK approach are quite reasonable. Apart from it is necessary to recognize that expressions for 
the Green functions are much more simple in comparison with any traditional regularization. 

In the framework of traditional regularizations the chiral-symmetry-breaking phase always co- 
exists with the chiral-symmetric phase, which is energetically unfavored. The existence of the 
symmetric phase is not evident for the implicit regularization in the original formulation of work 
[1]. But the transition to the symmetric phase with M = (in chiral limit mo = 0) can be easily 
performed with the proposed differential formulation of the implicit regularization. Really, taking 
M = in equation (8), we obtain = — and, therefore, l2{p'^) = —^^log^^ instead of 

(9). Then Ii = —-^^ at M ^ and from equations (12) and (13) it follows that /tt ^ and 
c ^ in correspondence with the trivial physics of the symmetric phase. 

An essential difference from other regularizations manifests the properties of Green functions 
in the Euclidean momentum region. In the BDK approach the meson propagators possess a non- 
physical singularity - a pole at the negative momentum square (Landau pole, or Landau ghost). ^ 
The meson propagators in the NJL model are scalar and pseudoscalar parts of the two-particle 
amplitude and can be written as 

" ~mo/M - 4z5nc(4M2 - p^)l2(p^) ^^^^ 

for the sigma-meson and 

n = !^ 

mo/M + Ugricp'^hip'^)' ^ ' 

for the pion. In the chiral limit the equation for the Landau pole, as it follows from equation (14), 
is 
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(20) 



where z = ^\ — At above values of the model parameters {ric = 3, fn = 93 MeV, M = 

468 MeV) this equation has solution p| = -4.29M^ = -(969 MeV)"^, i.e., Landau mass value 



^The existence of such pole was discovered firstly in quantum electrodynamics [7] and presents a characteristic 
feature of theories without an asymptotic freedom in deep-Euclidean region. The existence of the Landau ghost in a 
system of fermions coupled to a chiral field has been observed in work [8] 
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= —pj^ approximately twice larger in comparison with the quark mass: = 2.07M. Beyond 
the chiral limit (at mo / 0) this value changes very small (about 1%) due to the smallness of 
current quark mass niQ. 

The existence of the Landau pole is a serious problem of the BDK approach. In particular, any 
calculations with meson loops become problematic. Though such calculations exceed the frame- 
work of the one-loop approximation, which is a subject of work [1], the impracticability of these 
calculations means a principal impossibility of calculations of corrections to leading approximation 
and cannot be acceptable. 

Further, the Landau pole in quantum electrodynamics due to the smallness of fine structure 
constant a is located in the very distant asymptotic Euclidean region {{M1)^^^ ~ — mgexp{^}, 
where nie is the electron mass and a ~ 1/137), and its presence can be in principle ignored. Really, 
at the much smaller energies the quantum electrodynamics becomes a part of an asymptotically 
free grand unification theory with self-consistent asymptotic behavior. In the NJL model with 
the implicit regularization the Landau pole is located near the physical region of the model, and 
similar reasoning is impossible in principle. From the other hand, the traditional regularizations, 
such as the four-dimensional cutoff or the Pauli-Villars regularization, are free on this problem - 
the meson propagators in these regularizations do not have the Landau poles. Therefore, implicit 
regularization, possessing the certain appeal and simplicity, contains serious defect as the nearby 
Landau pole. 

To improve the situation a compromise approach is necessary, which conserves main features 
of implicit regularization and at the same time solves the problem of the Landau pole. Such 
compromise can be achieved with the Feynman regularization for logarithmically divergent integral 

^'^^'^ " / '{2^{ (M2 - (p + g)2)(M2 - g2) - (M2 -{p + qf){M^ - q') 1' ^^^^ 

where Mj. is a regulator mass (M^ > M^), and the definition of quadratically divergent integral I\ 
as before is made by relation (4). 

In Euclidean region < integral (21) can be represented as 

72 = ^F(-^) 
^ 167r2^ M2^' 

where 

N /"^ , M^/M"^ + uil-u)x , , 

F{x)= dulog^^ — ^- — 22 

^ ^ Jo l + u{l-u)x ^ ' 

Prove the absence of zeroes of this function at a; > 0. Evaluat ing eleme ntary integrals in equation 

(22) and introducing variables v = ^(-^1-1-1 — 1) and Vr = ^{^/^ + ^^~^) obtain the following 
expression: 

F = 2 [log -—^ + Vr log(l + —)-v log(l + -)1 . (23) 

1 + V Vr V 

If > M^, then Vr > v, and from elementary inequalities log > and Vr\og{l + ^) > 
'i;log(l -|- i) it follows that F > in the Euclidean region, i.e. the meson propagators do not 
possess the Landau pole with this definition. For other aspects this modified regularization is 
similar to implicit regularization. For the modified regularization 

^=(0) = -y^log^. (24) 

and, consequently, Ii is defined by same formula (11) with substitution 

Mo ^ Mr, 
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i.e. integration constant Mq everywhere is substituted by regulator mass M^. Such re-dcfinition 
of Ii still implies an existence of additional parameter M^, which is the integration constant of 
equation (4). 

Equation (15), which determines quark mass M, has exactly the same form for this modified 
regularization, and, consequently, all parameter values are the same. Remark, condition > 
is the automatical consequence of formula (12). Therefore, the proposed modification conserves 
main features of the implicit regularization and simultaneously solves the problem of Landau pole.^ 

A general receipt for definitions of divergent integrals (2) and (3) can be formulated as follows: 
integrals (3) are defined by formulae 

{h; l^; in = {h; I^; in'^'^'iM') - {h: I^; /fl^^^lM^). (25) 

Upper index BDK means that each integral with mass M or is defined by formulae of work 
[1].5 

Then integrals (2) are defined by equations (4) and (5). A consistency relation is the coincidence 
of integration constants, and, therefore. 

If = -jfc^/i (26) 

as in the BDK approach with consistency conditions. 

It should be noted, that due to linearity of consistency constraints considered in work [1], the 
analysis of symmetry conservation, which was made in this work, can be carried to the proposed 
modification. 

A principal distinction between the implicit regularization and the proposed ghostless modifi- 
cation consists in the treating of integral l2- Parameter Mq in the implicit regularization is the 
integration constant whereas parameter in the Feynman regularization (21) is introduced "by 
hands" and seems to be an outside object. The question arises, if really this object is some "secret 
parameter" , which ensures the self-consistency of the N JL model on the quantum level? We cannot 
answer to this question now with the full confidence, but the structure of the implicit regularization 
gives us an allusion to the possible realization of this idea. The term with regulator mass My. in the 
definition of integral I2 by formula (21) can be interpreted as the rise of some repulsive interaction 
(i.e. interaction with a negative coupling) at the one-loop level. The coupling constant of the NJL 
model is connected with the quark mass and the condensate by formula (17). If the quark conden- 
sate is negative the negative coupling corresponds to the negative quark mass. In this connection 
it is not out of place to recall about the negative solution of BDK equation (15) ("quark ghost"). 
In the physical case this solution is — 2M, where M is the dynamical quark mass. The existence 
of this solution can be an indication of such repulsive interaction, which becomes apparent on the 
one-loop quantum level. At given above parameter values we have \Mr\ 1.3M. This value, of 
course, is far from the expected value. When the absolute value of the quark condensate decreases 
ratio \Mr\/M arises, and at c = —210 MeV it equals to ~ 1.5. At c = —150 MeV it reaches 
expected value 2, but the last value of the condensate is phenomenologically unacceptable since it 
leads to the very large value of the current quark mass mo = 21 MeV. Therefore, the question 
of the possible "inter-termination of ghosts" in the framework of the NJL model remains to be 
open. 

* Another way to remove the Landau ghost in the NJL model is the method of work [9], based on the Kallcn 
Lehmann representation. This method has been applied in work [10] to the chiral cr-modeL However, the proposed 
receipt seems to be preferable in the framework of the implicit regularization since this regularization has deal with 
a set of integrals while the method of works [9, 10] should be applied to the proper two-point functions. 

®With such definitions integrals I2 and are convergent without additional regularization. Since for self- 
consistency of the receipt it is necessary to verify the coincidence of usual calculation of these integrals with the 
calculation by formula (25), based on the BDK definition of divergent integrals. Simple calculation demonstrates the 
reality of such coincidence. 
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In conclusion, discuss briefly the definition of model parameters beyond the chiral limit. At 
mo 7^ from eq.(19) it follows the pion spectrum equation 

mo = — m^-him^). (27) 

Practically this equation is used for definition of current quark mass mo- 

Taking approximation /2(m^) ~ -^2(0), we obtain the well-known current-algebraic Gell-Mann- 
Oakes-Renner formula 

mlfl ~ -mo < V^V' > • (28) 

This relation, of course, is regularization- independent (and what is more, model-independent). At 
rriTT = 135 MeV, c = —250 MeV we obtain from eq.(28) niQ = 5.2 MeV. The exploiting of 
regularization-dependent exact value Io{m^) changes this result very slightly (approximately of 1 
percent) due to the smallness of pion mass, and such specification, of course, lays out the framework 
of the model, it being right as for the implicit regularization as the modified variant. 

The sigma-mcson mass in the chiral limit is 2M. Beyond the chiral limit the equation on the 
sigma-meson mass follows from equation (18): 

mo = ^^^-{ml-AM').l2iml). (29) 
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With rather crude approximation l2{m'^) ~ -^2(0) and with taking into account equation (28) we 
obtain the well-known regularization-independcnt formula (see, e.g., [2]) 



mi c± 4M^ + m;. (30) 

The much more exact (but regularization-dependent) approximation is himl) ~ 72(4M2). In this 
case, we obtain 

ml ~ + Q.llml (31) 

for the implicit regularization, and 



ml ~ + 0.26m^ (32) 

for the modified variant. 

As we see from the above formulae, in the physical region of energies about 1 GeV the quanti- 
tative difference of these regularizations is not very significant. 
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